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1. INTRODUCTION 

If (M, g) is a Lorentz manifold and S a spacelike hyper surface, let h and k denote the induced 
metric and second fundamental form of S respectively. The mean curvature of S is the trace tr^A;. 
An interesting class of spacelike hypersurfaces are those whose mean curvature is constant (CMC 
hypersurfaces). The Lorentz manifolds of primary interest in the following are those which possess a 
compact Cauchy hypersurface and whose Ricci tensor satisfies r{V, V) > for any timelike vector V . 
They will be referred to as cosmological spacetimes. The curvature condition is known as the strong 
energy condition and implies certain uniqueness statements for CMC hypersurfaces. Under very 
general conditions, if a region of a cosmological spacetime is foliated by compact CMC hypersurfaces, 
then each of these has a different value of the mean curvature and a time coordinate t can be defined 
by the condition that its value at a point be equal to the mean curvature of the leaf of the foliation 
passing through that point. A time coordinate of this type will be referred to in the following as a 
CMC time coordinate. 

The questions, whether a cosmological spacetime contains a compact CMC hypersurface and 
which values the mean curvature can take on a hypersurface of this kind may seem purely geometrical 
in nature. However it turns out that the answers to these questions depend crucially on factors which 
have no obvious geometrical interpretation, but which have a physical meaning, when the Lorentz 
manifold is considered as a model for spacetime. 

The Einstein equations for a Lorentz metric g take the form G = 8ttT, where G is the Einstein 
tensor of the metric g and T is the energy-momentum tensor. To get a determined system of evolution 
equations for the geometry and the matter, it is necessary to say more about the nature of the matter 
model used. This means specifying some matter fields, denoted collectively by 0, a definition of T in 
terms of g and <p an d the differential equations which describe the dynamics of the matter. Putting 
these things together with the Einstein equations gives a system of evolution equations, the Einstein- 
matter equations. It will be seen that the existence of global foliations by CMC hypersurfaces in a 
cosmological spacetime which is a solution of the Einstein-matter equations depends essentially on 
the matter model chosen. 

Let s denote the scalar curvature of the metric g. The spacetime is said to satisfy the weak energy 
condition if r(V,V) > (l/2)sg(V, V) for all timelike vectors V. (Note that, despite the terminology, 
the strong energy condition does not imply the weak one.) There is a topological obstruction to the 
existence of a compact spacelike hypersurface with vanishing mean curvature (maximal hypersurface) 
in a spacetime which satisfies the weak energy condition. In a cosmological spacetime this implies 
that if the Cauchy hypersurface is a manifold which admits no Riemannian metric with non-negative 
scalar curvature then the spacetime contains no compact maximal hypersurfaces. Moreover, there are 
strong restrictions in the case of a manifold which admits no Riemannian metric of positive scalar 



curvature. The main questions are whether a cosmological spacetime satisfying the weak energy 
condition and containing at least one compact CMC hypersurface can be covered by a foliation 
of compact CMC hypersurfaces and whether the mean curvature of these hypersurfaces takes on 
all values not forbidden by the topological obstruction already mentioned. In other words, does it 
take all values in the interval (— oo,0) or (0, oo) in the case where there is an obstruction and all 
real values in the case there is none? This statement concerning the values attained by the mean 
curvature is equivalent to a statement of whether a solution of the Einstein-matter equations exists 
globally in a CMC time coordinate. The theorem proved in Section 2 shows by example that the 
answer is negative if no restriction is put on the matter model used. In Section 3 results are reviewed 
which show that under certain symmetry assumptions there are matter models for which the answer 
is positive. Possible extensions of these results are also discussed. For general information on the 
points just mentioned, the reader is referred to [1], which is complementary to the treatment here. 

These results are related to the concept of crushing singularities. A cosmological spacetime 
is said to have an initial crushing singularity if there is a compact spacelike hypersurface S to and 
a foliation of the past of S to by compact spacelike hypersurfaces S t , such that as t tends to its 
limiting value towards the past, the mean curvature of St tends uniformly to — oo. (There is a similar 
definition with 'initial' replaced by 'final', 'past' replaced by 'future' and — oo replaced by oo.) If 
the S t form a CMC foliation and the parameter t labelling the leaves of the foliation is a CMC time 
coordinate, then a sufficient condition in order to have an initial crushing singularity is that t should 
take all values in the interval (— oo, to). The positive results of Section 3 prove that this is true 
under certain hypotheses on the matter model and the symmetry of the spacetimes considered, and 
so prove the existence of crushing singularities under certain circumstances. On the other hand, the 
results of Section 2 provide examples of spacetimes where an initial (or final) crushing singularity is 
not present. Informally, this can be expressed by saying that the initial (or final) singularity in these 
spacetimes is not crushing. It will now be indicated briefly how the statement about the absence of 
crushing singularities follows from the theorem proved in Section 2. If there is an initial crushing 
singularity, the hypersurfaces St provide barriers, as they are used in the well known existence theory 
for CMC hypersurfaces. Thus it can be assumed without loss of generality that the S t are in fact CMC 
hypersurfaces and that t is a CMC time coordinate. Combining this with the uniqueness theorems 
for CMC hypersurfaces would show that the local CMC foliation which exists close to the initial 
hypersurface could be extended to arbitrarily negative values of t, contradicting the conclusions of 
the theorem. 

The techniques which are used to obtain the positive and negative results are closely related. 
They will be explained in the case of spacetimes with U(l) x U(l) symmetry where the Cauchy 
hypersurface has the topology of a three-dimensional torus. Assume that a solution of the Einstein- 
matter equations satisfying the strong energy condition with a CMC Cauchy hypersurface of topology 
T 3 is invariant under the action of the group U(l) x U(l) consisting in rotating two of the three S 1 
factors. Then (see [2]) in a neighbourhood of the initial hypersurface the metric can be written in 
the form: 

-a 2 dt 2 + A 2 [(dx + f3 l dt) 2 + a 2 ~g AB (dy A + (3 A dt)(dy B + (3 B dt)\ (1.1) 

where t is a CMC time coordinate. Here the coordinates are t,x,y 2 ,y 3 . Upper case Roman indices 
take the values 2,3 while lower case ones take the values 1,2,3, the value 1 corresponding to the 
coordinate x. The functions a, f3 a , A and g AB depend on t and x and g AB has unit determinant. 
They are periodic in x. The quantity a depends only on t. Some of the field equations are: 

d 2 x (A l/2 ) = -i^ 5/2 [f(ifi - \tf - It 2 + 2 VA r, A + k AB R AB + X AB \ AB + 16tt P ] (1.2) 
d 2 a + A~ 1 d x Ad x a = aA 2 ^^ - ±t) 2 + \t 2 

+ 2r] A r] A + k AB k AB + 4?r(p + trS)] - A 2 (1.3) 

d x K x + 3A~ 1 d x AK 1 - A~ 1 d x At - R AB \ AB = 8ttJA (1.4) 
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8 t a = a[-8 x ^ + |a(3ifi - t)} 
d t A = -aK 1 A + d x ((3 1 A) 



(1.5) 
(1.6) 



The quantity K\ appearing in these equations is an eigenvalue of the second fundamental form. 
Alternatively it may be thought of as an auxiliary quantity defined in terms of the basic quantities 
contained in (1.1) by (1.5). The quantity 77^4 is given by 

7M = {l/2) a - 1 Aa 2 ~g A Bd x ^ B (1.7) 

while Rab and Xab are the tracefree parts of the second fundamental forms of the group orbits 
in spacetime corresponding to the normal vector to the hypersurface t=const. and the normal to 
the orbit in the hypersurface i=const., respectively. The components c/ab can be parametrized as 
follows: 

922 = e w coshV, 533 = e~ w cosh V, <?23 = sinhF (1.8) 
In terms of W and V the squares of Rab and Xab have the explicit forms: 

\ AB X AB = \ A~ 2 (cosh 2 VW 2 + V 2 ) (1.9) 
~k ab ~k ab = la" 2 [cosh 2 V(W t - P l W x ) 2 + (V t - pV x )*\ (1.10) 

The quantities p, J and tiS denote the energy density, the matter current and three times the mean 
pressure, respectively. 

In the proofs of the existence and non-existence theorems it is important to have estimates for a 
solution of these equations on a finite time interval (ti,t 2 ) with t 2 < in terms of the data it induces 
at some intermediate time to- Consider a point where a attains its maximum on a hypersurface of 
constant time. The strong energy condition implies that p + trS > and so it follows from equation 
(1.3) that a < 3/t 2 . In fact this estimate is true, and can be proved in the same way, without any 
symmetry assumption on the spacetime. The next step is to use a generalization of an argument of 
Malec and O Murchadha [3] to show that if the dominant energy condition holds then \K\\ < 5\ti\, 
l^ -2 ^^!! < 2\t\\. (For details of this argument see [2].) With these estimates (which depend very 
much on the symmetry assumption) in hand, equations (1.5) and (1.6) can be integrated in time to 
give the following bounds for a and A. 

a(t ) exp(-Cft - t \) < a(t) < o(t ) exp(C|t - t \) 

sup{A(t,x),A-\t,x)} < sup{p(t )|| oo ,p- 1 (t )|| oo }exp(exp(C|t-t |)) 

Along the way it also comes out that l^/? 1 ! satisfies the same kind of upper bound as a. For details 
the reader is once again referred to [2]. These bounds are probably far from optimal but they are 
sufficient for the present purposes. The general point is that many aspects of the geometry can 
be bounded in terms of t±, t 2 and the maximum and minimum values of a and A on the initial 
hypersurface t = to- Using this information and integrating equation (1.2) in space at each fixed 
time gives a bound for the integral J Q n Q(t, x)dx where 

Q = 167rp + 2r] A r] A + k ab R A b + \ AB \ab (1.12) 

which is uniform in time on the given interval. 

2. EXAMPLES OF SINGULARITY FORMATION 

A matter model which is well known for its tendency to form singularities is dust. In this section 
a type of example will be presented which shows that global existence in a CMC time coordinate 
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does not hold for the Einstein-dust system. The main interest of this is the contrast it provides with 
the known positive results for some other matter models which will be reviewed in the next section. 
To put it another way, this example says something about the sharpness of the results previously 
obtained. The intuitive idea behind the example is very simple. If two dust particles start out close 
together and with velocities which are moderate in magnitude and opposite in direction then they 
should collide after finite time producing a shell-crossing singularity. To turn this intuition into a 
proof we can use the bounds on various geometric quantities presented in the last section. However 
these estimates alone do not suffice. The difficulty is connected with the phrase 'in finite time'. The 
problem is that a finite amount of proper time along the wordline of a particle might correspond to 
an infinite amount of coordinate time. This would correspond to the phenomenon known as 'collapse 
of the lapse', where the lapse function a becomes very small. To show that this does not happen it 
is necessary to have a positive lower bound on the lapse function on a finite interval of CMC time. 
This is the subject of the following lemma. 

LEMMA (supporting the lapse) Consider a solution of the Einstein- matter equations with £7(l)x[/(l) 
symmetry on a finite interval (ii, £2) of CMC time with t 2 < and let to belong to this time interval. 
Suppose that the strong and dominant energy conditions hold. Then there is constant C, depending 
only on t\, t 2 and the maximum values of a, a -1 , A and A -1 on the hypersurface t = t such that 
a" 1 < C. 

Proof. Equation (1.3) can be written in the form 

d x (Ad x a) = aA 3 P - A 3 (2.1) 

where P is non-negative. Under the given hypotheses G = J„ 27r P(x)dx is bounded. Consider now a 
fixed time t and an interval I = [x\, x 2 ] of values of x. Let L = \x 2 — x\\ and denote by F and F' 
the maximum values of Aa and d x (Aa) respectively on the interval /. Now for any x G I 

(Aa)(x) = (Aa)( Xl ) + f d x (Aa)(y)dy 

J Xl 

and the integral can be bounded by LF'. Hence: 

F < (Aa)(xi) + LF' (2.2) 

On the other hand 

d x (Aa)(x) = (Ad x a)(x) + {A~ x d x A{Aa)){x) 

= (Ad x a)( Xl ) + f d x (Ad x a)(y)dy + (A" 1 ^A(^a))(x) 

= d x (Aa)( Xl ) - (A^d^Aa))^) + f {aA 3 P - A 3 )(y)dy + {A~ 1 d x A(Aa)){x) 

J X\ 

In the last line equation (2.1) has been used. Since A~ l d x A has already been bounded, the second 
and fourth terms in the final expression are bounded by CF for some positive constant C. Consider 
now the third term. The second term in the integrand makes a negative contribution, while the first 
can be bounded by CFP. Hence: 

F' < d x (Aa)(x\) + C(l + G)F (2.3) 
It follows by substituting (2.3) into (2.2) that 

F < (Aa)( Xl ) + L[(d x {Aa){ Xl ) + C(l + G)F] (2.4) 
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If L < \[C{1 + G)]~ l then this implies that 

F < 2[(Aa)( Xl ) + Ld x (Aa)( Xl )} (2.5) 
Similarly, under the same restriction on L: 

F' < 2[d x (Aa)( Xl ) + C(l + G)(Aa)( Xl )} (2.6) 

Thus we have the inequalities 

(Aa)( X + L) < 2[{Aa){ X ) + Ld x (Aa)( X )} 
d x {Aa)(x + L)< 2[d x {Aa){ X ) + C(l + G)(A*)(x)] ' 

Now two cases will be considered. The first is where C(l + G) < (47r) _1 . Then L can be chosen to 
be 2n and it follows that if Aa takes its minimum at a given time at x_ then 

{Aa){x) <2{Aa){x-) (2.8) 

If C(l + G) > (47r) _1 then it is possible to divide the circle into a number of equal intervals, starting 
at X - whose length L satisfies the desired inequality and whose number does not exceed C(l + G). 
Let z( X ) = max{(Aa)(x),d x (Aa)(x)}. Then 

z(x + L) < 2C(1 + G)z(x) (2.9) 

Applying this on k successive intervals gives 

z{ X + kL) < 2 k C k {l + G) k z( X ) (2.10) 

On the other hand Aa( X ) can be bounded by a fixed constant times z( X + kL) for some k < C(l + G). 
There results an estimate of the form 

(Aa)(x) < [C(l + G)] c{1+G) (Aa)( X .) (2.11) 

As discussed in the previous section, the assumptions of the lemma imply bounds on A and A -1 . 
Hence (2.8) and (2.11) imply an estimate of the form || 

c^lloo Co_, where o_ is the minimum value 
of a at the given time. On the other hand, integrating equation (2.1) with repect to x shows that: 

»2tt i-2-it 



/*Z7T />Z7T 

/ A 3 (x)dx = / {aA 3 P)(x)dx < C\\a\ 
Jo Jo 



Since f Q ^ A 3 ( X )d X is bounded from below, putting these estimates together gives the desired lower 
bound for a_ . 

The estimates (2.8) and (2.11) resemble the specialization of Harnack's inequality to one space 
dimension with the difference that the constant in the inequality only depends on the L 1 norm of P 
rather than its L°° norm. 

Next the above lemma will be applied to the Einstein-dust system. The matter fields are a non- 
negative function fi, the proper energy density of the dust, and a unit vector u a , the four- velocity 
of the dust particles. The energy-momentum tensor is given by T a p = \iu a up and the equations 
describing the dynamics of the matter fields are simply V Q T a ' 3 = 0, which is of course a necessary 
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compatibility condition for the Einstein equations. The dominant, strong and weak energy conditions 
are satined by this matter model. The theorem to be proved is the following: 

THEOREM Let to be a negative real number and e > 0. Then there exist initial data with constant 
mean curvature to for the Einstein-dust system such that in any Cauchy development of these data 
the CMC foliation which exists in a neighbourhood of the initial hypersurface cannot be extended to 
values of the mean curvature greater than to + e- Similarly there exist data for which the foliation 
cannot be extended to values of the mean curvature less than to — e. 

Proof. The initial data to be constructed will have U{\) x U(l) symmetry and so it is possible to 
use the form (1.1) of the metric. They will also be such that the initial velocity of the dust particles 
is in the x-direction. Then the velocity can be parametrized by the inner product v of the four- 
velocity u a with the unit vector A~ 1 d/dx. The integral curves of u a are geodesies and hence there 
are conservation laws corresponding to the Killing vectors d/dy A . This means that if the velocity is 
initially in the x-direction, it remains so. Thus these dust solutions can be described completely by 
the two functions fj, and v. 

Initial data for dust spacetimes can be constructed using the conformal method. (A general 
discussion of this method can be found in [4].) However the situation here is sufficiently simple that 
we will not need all of the general theory. The data constructed will be such that Rab = 0, Xab = 
and t\a = 0. The matter quantities occurring in the constraints are related to the matter fields fi 
and v in the special case of the type of data under consideration by: 

P = Kl + v 2 ) 
J = fi (l + v 2 ) 1 / 2 v 

To construct initial data, choose first the constant value t = to for the mean curvature and a non- 
negative function jl, a function v and a scalar function K\(x) on the circle. The solution of the 
constraints is sought in the form: 

K x = \t + A~*(K x - \t) 
H = A~ 4 fi 

and the Lichnerowicz equation, which is just (1.2) rewritten in terms of the rescaled quantities, reads: 
d 2 x {A 1 / 2 ) = -^[lA-^K, - \tf - It 2 + 167rA- 4 /i(l + v 2 )] 

In terms of the rescaled quantities the equation (1.4) takes the simple form d x (K 1 — |i) = 8ttJ, 
where J = p,(l + This equation can be solved if and only if J has integral zero. One way of 

ensuring that this condition is satisfied is to impose the symmetry conditions that fl(x) = /2(7r — x) 
and v(x) = —v(tt — x). When this equation is solvable the L°° norm of K\ — |i can be estimated in 
terms of the L 1 norm of J. The next step is to solve the Lichnerowicz equation and for that an idea 
will be borrowed from the general method, namely that of using sub- and supersolutions. In order 
for this to run as smoothly as possible, assume that fi is bounded below by a positive constant B. 
Then it is possible to find constant sub- and supersolutions, namely 

A_ =to 1 (24vr J B) 1/2 

(2 14) 

A + = max{t - 1 (48vr||/i|| oo (l + IMlL)) 172 ^!^ " ^llL) 2/3 } 

These ensure the existence of a solution of the Lichnerowicz equation and give pointwise estimates for 
the solution. This allows A and A -1 to be bounded pointwise in terms of the L°° norms of K\ — 
\x and v and the constants to and B. 
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(2.12) 



(2.13) 



The above provides a possibility of constructing a variety of initial data in such a way that the 
quantities entering the hypotheses of the lemma above can be easily controlled in terms of the free 
data. To obtain the desired example a specific subclass of this data will be considered. Let x\ < x 2 
and choose initial data for v such that v\ = v{x\) = 1 and V2 = v{x2) = — 1. Consider two dust 
particles which start at time to at the points x\ and X2 and with the velocities v\ and V2 respectively. 
Let the positions and velocities of these particles at time t be denoted by xi(t), X2(t), vi(t) and V2(i). 
These quantities satisfy the equations: 

dx/dt = aA^v/y/l + v 2 - (3 1 (2.15) 
dv/dt = -A _1 aVl + t; 2 + aKv (2.16) 

It is elementary to see that on a given finite interval of CMC time bounds for v\ and V2 can be obtained 
using GronwalPs inequality. Moreover these bounds do not depend on the distance |a?2(io) — £i(*o)|- 
The idea now is to assume that for a family of initial data of this type with \x2(to) — #i(io)| tending 
to zero the corresponding solutions exist at least up to and including the time to + e and to show 
that this assumption leads to a contradiction. Note that if to + e > the statement of the theorem 
is an immediate consequence of the non-existence of maximal hypersurfaces. Thus it is assumed in 
the following that to < — e. 

Let C be a constant greater than 5sup to<t<to _|_ e {(||A _1 a / || 00 + HaKjloo)}. By what has been 
said above we know that the data can be chosen so that a single constant C works for all data in 
the family. It will now be shown that for the solution evolving from any one of these data v\ (t) > | 
for all times t < to + e such that to < t < to + C~ l . Let t* be the largest time in the interval 
[to, to + e] such that Vi(t) > ^ for all t in the interval [<o,i*]- If i* < h + C~ x and t < t + e, let 
t' be the last time before i* that v\{t) was equal to unity. Now v(t*) = \. On the other hand it 
follows by integrating (2.16) from t' to t* that v(t*) > |. This contradiction shows that in fact either 
t* > t + C _1 or t* = to + e. This gives the desired conclusion. In a similar way it can be shown 
that v 2 (t) < — \ for all times t < t + e such that t < t + C~ x . For convenience of notation, let 
t 3 = min{t + e, t + C" 1 }. From (2.15): 

d/dt(xi - x 2 ) = aA~ l v 1 /^ 1 + v\ - aA~ 1 v 2 /^ 1 + 4 - P 1 (x 1 ) + p 1 (x 2 ) (2.17) 

On the interval [0, t 3 ] we h ave a lower bound for v\, an upper bound for f 2 and a crude upper bound 
for ^\ + v\ and \J\ + v\. Moreover, we have a lower bound for aA~ l . (It is at this point that the 
lemma on supporting the lapse is used.) Putting all this together gives a negative upper bound for 
the sum of the first and second terms on the right hand side of (2.17). On the other hand 

\(3\x 1 )-l3\x2)\ < R/^Uooki-^l (2.18) 

Hence by choosing \x\{to) — x 2 {to)\ small it can be ensured that the sum of third and fourth terms on 
the right hand side of (2.17) is negligeable in comparison with the sum of the first and second terms. 
Hence it can be arranged that x\(t) — x 2 (t) tends to zero after a time smaller than t 3 . However this 
contradicts the existence of a regular solution of the Einstein-dust equations up to time to + e, since 
in a regular solution the world lines of dust particles can never cross. Thus the first statement of the 
theorem has been proved. The proof of the second statement is strictly analogous. 

3. EXAMPLES OF GLOBAL REGULARITY 

There are matter models for which the situation is quite different from that presented in the last 
section, in that for spacetimes with U(l) x U(l) symmetry possessing a CMC Cauchy surface global 
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existence in CMC time is obtained. In this context 'global existence' means existence on the longest 
time interval consistent with the topological obstruction discussed in Section 1. A global existence 
theorem in this sense was proved in [2] for collisionless matter described by the Vlasov equation and 
for the massless scalar field (or more generally wave maps) . The nature of these two matter models 
will now be recalled. In the case of collisionless matter, the matter field is a non-negative real-valued 
function / on the space of future-pointing unit timelike vectors in spacetime (the mass shell) and 
the equation which describes the dynamics of the matter is the Vlasov equation. This simply says 
that the function / is constant along the curves which are the natural lifts of timelike geodesies to 
the mass shell. The energy-momentum tensor at a given spacetime point is obtained by integrating 
the product of / with a suitable weight over the part of the mass shell over that point. Details can 
be found in [5]. In the case of the massless scalar field, the matter field is a real- valued function 
<fi on spacetime which is supposed to satisfy the wave equation. The energy-momentum tensor is 
of the form T = d(f> (g> d(f) — ^\d(j)\ 2 g. A wave map, which is a Lorentzian analogue of the harmonic 
maps familiar in Riemannian geometry, is a generalization of this, where the field, instead of taking 
values in the real numbers, takes values in an arbitrary complete Riemannian manifold, known as 
the target manifold. It is worth to note that the vacuum case, i. e. the case T = 0, is contained 
in these results as the special case / = or <\> = 0. The global existence problem is hard even in 
the case of U(l) x U(l) symmetric vacuum spacetimes. The physics issue is whether gravitational 
waves propagate smoothly for arbitrarily long times or whether they could develop shocks. It is 
reasonable to expect that when a global existence theorem in CMC time holds, it should not be 
possible to extend the spacetime beyond the region covered by the CMC foliation, while maintaining 
the property that the initial hypersurface be a Cauchy hypersurface for the extended spacetime. 
Unfortunately, this statement has not been proved up to now. 

The example of collisionless matter is particularly interesting due to the fact that dust solutions 
can be considered as distributional solutions of the Einstein- Vlasov system [6]. They have Dirac 5- 
function dependence on the velocity variables. The results already quoted show that approximating 
the 5-function by smooth functions in the initial data leads to a dramatic change in the long-time 
behaviour. 

The basis of the global existence theorem is formed by the estimates mentioned in Section 1. 
One then proceeds by bounding higher and higher derivatives of the metric and the matter fields on 
the given finite interval (ti,t 2 ). When all derivatives have been bounded it follows that the solution 
can be extended to the closed interval. A local existence theorem (which is a consequence of standard 
results on the Cauchy problem and the existence of CMC foliations) then allows it to be extended to 
a longer time interval. Finally, consideration of the maximal interval of existence implies the desired 
global theorem. 

The fact that higher derivatives can be bounded is connected to the fact that the given matter 
model does not form singularities in a given regular spacetime. This condition is violated by dust. 
The quantity p(t, x)dx is a bounded function of t on the given time interval, independent of the 
matter model. On the other hand, it is to be expected that p can blow up pointwise on a finite 
time interval in the case of dust. For other matter models it could happen that, although the energy 
density remains bounded, regularity breaks down at a higher level. This could result from formation 
of shocks by the matter. Coming back to the general case, the PDE problem to be studied is that of 
a semilinear equation of wave map type for W and V coupled to the matter equations. The target 
space of the wave map is the hyperbolic plane. This system is defined on a curved two-dimensional 
geometry defined by a, (3 1 , a and A. What needs to be shown is that the finiteness of a certain 
C k norm of the two-dimensional geometry implies that of a similar norm of W, V and the matter 
quantities. Conversely, the finiteness of the latter implies, via equations (1.2)-(1.6), the finiteness of 
a stronger norm of the two-dimensional geometry. This allows higher derivatives of all quantities of 
interest to be bounded inductively in favourable cases. The details of the argument for bounding low 
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order derivatives depend very much on the particular matter model. 

There are similar results for spherically symmetric spacetimes with a Cauchy hypersurface of 
topology S 2 x S 1 . In that case there is no obstruction to the existence of a maximal hypersurface 
and all real values are attained by the mean curvature [7,8]. Moreover, it can be shown that the 
CMC foliation covers the entire spacetime. There are also some other cases where results are known 
but they do not involve larger classes of solutions than those discussed above (for details see [1]). 
All the known positive results are for spacetimes with at least two symmetries, so that the problem 
effectively reduces to studying a system in one (or less) space dimension. It seems that even the 
case with one symmetry is very difficult, not to mention the general case. However, investigations 
into generalizing the results reported here to those cases are being carried out. It should be noted 
that symmetries with fixed points also lead to difficulties so that, for instance, the case of spherical 
symmetry on S 3 remains open. 

Another possible direction for generalizations is to keep the high symmetry but to relax the 
assumptions on the matter model. There are two basic types of matter model to be considered. 
There are the phenomenological matter models (e.g. dust, collisionless matter) and the field-theoretic 
matter models (e.g. massless scalar field, wave maps). The phenomenological matter models represent 
a macroscopic description of matter. They often have a tendency to form singularities in finite time 
in a given smooth spacetime. Dust is a good example. It would probably be possible to prove an 
analogue of the theorem of Section 2 for perfect fluids with pressure. The difference would be that, 
since the expected singularities are shocks, the energy density would probably remain finite at the 
time when the CMC foliation broke down. To try and get positive results for a fluid, it would be 
natural to introduce viscosity. Unfortunately, the concept of viscous fluids is known to be problematic 
even in special relativity. A case where it is difficult to make predictions is that of the Boltzmann 
equation. This is because, to the author's knowledge, there is up to now not even a global existence 
theorem for classical solutions of the special relativistic Boltzmann equation in the case of spatially 
homogeneous initial data. 

Field theoretic matter models are supposed to describe matter at a more fundamental level. 
They seem less prone to forming singularities in a given smooth spacetime than phenomenological 
matter models. They do have a different kind of problem, which may purely be an incompatibility 
with the known techniques, rather than an essential difficulty. The problem is that in certain steps 
of the proofs one would like to have the non-negative pressures condition, i.e. the condition that 
T(X,X) > for all spacelike vectors X. This condition is almost never satisfied by field-theoretic 
matter models. It is not even satisfied by the massless scalar field. However in the case of[/(l)x[/(l) 
symmetry (or spherical symmetry) it suffices that T(X, X) > for spacelike vectors X orthogonal 
to the orbits of the symmetry group. This latter condition is satisfied by the massless scalar field 
and, more generally, by wave maps. It is not satisfied by an electromagnetic field, a Yang-Mills field 
or a massive scalar field. This also has the inconvenience that the positive results on existence of 
CMC hypersurfaces in the case of collisionless matter do not obviously extend to the case of charged 
collisionless matter, coupled to an electromagnetic field. An argument used in [2], and which does 
not require non-negative pressures might allow one to circumvent that difficulty in the U(l) x U{1) 
symmetric case. However, it does not, as it stands, apply to the spherically symmetric case. The 
massive scalar field is even worse, since it does not satisfy the strong energy condition, which is the 
standard condition used to ensure uniqueness of CMC hypersurfaces. Perhaps some use can be made 
of the fact that the condition is satisfied in an average sense ([9], p. 95). It would be desirable to 
have more flexible techniques and clearly a lot remains to be learned in this area. 

To conclude, it is in order to make some remarks which put the question of the existence of global 
foliations by CMC hypersurfaces into a wider context. One of the most important mathematical 
open questions in general relativity is the cosmic censorship hypothesis of Penrose. It is closely 
related to the issue of the global behaviour of solutions of the Einstein equations corresponding to 
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general initial data. In [10] Eardley and Moncrief suggested that CMC hyper surf aces could be useful 
in trying to confirm this hypothesis. More generally, results of the kind discussed in this paper 
represent knowledge about long-time existence of solutions of the Einstein equations. This is not by 
itself enough to say something about cosmic censorship. For that one would not only need to know 
long time existence but also need detailed information on the asymptotic behaviour of the solutions. 
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